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We search for the factorization-suppressed decays B XcoK^*' and B — > Xc2K'*' , with XcO and 
Xc2 decaying into J/tpJ, using a sample of 124 x 10^ BB events collected with the BABAR detector 
at the PEP-II storage ring of the Stanford Linear Accelerator Center. We find no significant signal 
and set upper bounds for the branching fractions. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Nonleptonic decays of heavy mesons are not easily 
described because the process involves quarks whose 
hadronization is not yet well understood. The factoriza- 
tion hypothesis allows one to make some predictions ij 
by assuming that a weak decay matrix element can be de- 
scribed as the product of two independent hadronic cur- 
rents. Under the factorization hypothesis, B ccK^*'> 
decays are allowed when the cc pair hadronizes to J/ip , 
ipi^iS) or Xci, but suppressed when the cc pair hadronizes 
to XcO or Xc2 0]' Here, if^*^ represents either K or K* . 
In lowest-order Heavy Quark Effective Theory, there is 
no J > 2 current to create the tensor Xc2 from the vac- 
uum. The decay rate to the scalar XcO is zero due to 
charge conjugation invariance Q. 

Belle has recently observed XcoK^ decays with 

a branching fraction (BF) of ± 1.1) x 10"'' |1 

using XcO decays to tt+tt^ or K^K^ . BABAR has con- 
firmed the observation using the same decays with a 
branching fraction of (2.7 ± 0.7) x 10"^ 0, somewhat 
lower than, but compatible with, the Belle measurement. 
These results are of the same order of magnitude as the 
BF of the decay B^ XaK^ and are surprisingly 
large given the expectation from factorization. Using the 
hadronic XcO decays, CLEO has obtained an upper limit 
on XcoK° of 5.0 x 10"'' 0. Non- f actor iz able con- 

tributions to B^ XcoK^ decays due to rescattering of 
intermediate charm states have been considered theoret- 
ically 0, and similar branching fractions are predicted 
for decays to XcO and Xc2- No predictions are available 
for B decays to Xc(o,2)K*, but the branching fraction 
of decays to K* may be expected to be similar to the 
branching fraction of decays to K. The measurement 
of _B — > Xc(o,2)^*'*'' should improve our understanding of 
the limitations of factorization and of models that violate 
factorization. 

In this Letter we report a search for the decays 
B — > XcjK'^*\ J = 0,2, using the radiative decays 
XcJ J/i^l, with branching fractions of (1.18 ± 0.14)%, 
(20.2±1.7)%, respectively [g- Since the radiative branch- 
ing fraction for the XcO decay (including subsequent J/^p 
decay to £~^£~) is much smaller than the correspond- 
ing TT+vr^ or K'^K^ branching fractions, the search for 
the J5+ XcqK^ decay is less sensitive than previous 
searches, but it is free from the interference with the 
non-resonant decays to three mesons that affect the lat- 
ter. The data used in this analysis were obtained with 
the BABAR detector at the PEP-II storage ring, compris- 
ing an integrated luminosity of 112fb~^ of data taken at 
the T(4S') resonance. 



The BABAR detector is described elsewhere |9j. Sur- 
rounding the interaction point, a five-layer double-sided 
silicon vertex tracker (SVT) provides precise reconstruc- 
tion of track angles and i?-decay vertices. A 40-layer drift 
chamber (DCH) provides measurements of the transverse 
momenta of charged particles. An internally reflecting 
ring-imaging Cherenkov detector (DIRC) is used for par- 
ticle identification (PID). A CsI(Tl) crystal electromag- 
netic calorimeter (EMC) detects photons and electrons. 
The calorimeter is surrounded by a solenoidal magnet 
providing a 1.5-T field. The flux return is instrumented 
with resistive plate chambers used for muon and neutral- 
hadron identification. 

The channels considered here are B -> XcK^*^ with 
Xc J/i^l and J/'0 where £ is e or fi; K is A'+ 

or K° (^ TT+TT-); K*° K+n' or K°tt°; K*+ K+tt° 
or K'^TT^ ; and tt*^ —f 77. Charge-conjugate modes are in- 
cluded implicitly throughout this paper. Event selection 
is optimized by maximizing e/^/B, where e is the signal 
efficiency after all selection requirements and B the num- 
ber of background events, estimated with T{4S) — > BB 
and e+e~ — > qq Monte Carlo (MC) samples. 

Candidate J/'0 mesons are reconstructed from a pair 
of oppositely charged lepton candidates that form a good 
vertex. Muon (electron) candidates are identified with a 
neural-network (cut-based) selector and loose selection 
criteria. Electromagnetic depositions in the calorimeter 
in the polar-angle range 0.410 < 6iab < 2.409 rad that are 
not associated with charged tracks, have an energy larger 
than 30 MeV, and a shower shape consistent with a pho- 
ton are taken as photon candidates. For J/ip e+e~ de- 
cays, electron candidates are combined with nearby pho- 
ton candidates in order to recover some of the energy lost 
through bremsstrahlung. The lepton-pair invariant mass 
must be in the range [2.95, 3.18] GeV/c^ for both lepton 
fiavors. The small remaining background is mainly due 
to J/ip mesons not originating from Xc decays. 

We form candidates from oppositely-charged tracks 
originating from a common vertex with invariant mass 
in the range [487, 510] MeV/c^ The K° flight length 
must be greater than 1 mm, and its direction in the plane 
perpendicular to the beam line must be within 0.2 rad of 
the K'^ momentum vector. Charged kaon candidates are 
identified with a likelihood selector, based on information 
from the DIRC, and dE/dx in the SVT and in the DCH. 

A 7r° candidate is formed from a pair of photon can- 
didates with invariant mass in the interval [117, 152] 
MeV/c^ and momentum greater than 350MeV/c. K* 
candidates are formed from Kt: combinations with an 
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invariant mass in the range [0.85, 0.94] GeV/c^. 

The J 14': K^, and tt" candidates are constrained to 
their corresponding nominal masses to improve the 
resolution of the measurement of the four-momentum of 
their parent B-candidate. The Xc candidates are formed 
from J/ij) and photon candidates. The photon is required 
to have an energy greater than 0.15 GeV and not to be 
part of 7r° candidates in the mass range [0.125, 0.140] 
GeV/c2. 

Candidate B mesons are formed from Xc and if^*^ 
candidates. Two kinematic variables are used to fur- 
ther remove incorrectly reconstructed B candidates. The 
first is the difference /S.E = Eg — E^^^^^ between the 
-B-candidate energy and the beam energy in the T{AS) 
rest frame. In the absence of experimental effects, re- 
constructed signal candidates have Ai? = 0. The typical 
AE' resolution is 20 MeV for channels with only charged 
tracks in the final state, and 25 MeV, with a low /S.E 
tail due to energy leakage in the calorimeter, for chan- 
nels with a tt". The second variable is the beam-energy- 
substituted mass Toes = (^beam~-PB )^^^' whcrc p|j is the 
momentum of the B-candidate in the T(AS) rest frame. 
The energy substituted mass toes should peak at the B 
meson mass, 5.279 GeV/c^. Typical resolution for Ai? is 
2.7 MeV/c^. For the signal region, Ai? is required to be 
in the range [—35, -1-20] MeV for channels involving a 7r°, 
and within ±20 MeV otherwise. We require toes to be in 
the range [5.274, 5.284] GeV/c^. If more than one B can- 
didate is found in an event, the one having the smallest 
|Ai?| is retained. 

The observation of Xc2 could be complicated by the 
presence of the prominent Xci peak. This is mitigated by 
measuring the spectrum in the variable mg+i-^ — mi+g- . 
The efficiencies obtained from fits to the mass difference 
distribution for exclusive MC samples, where one B de- 
cays to the final state under consideration and the other 
inclusively, are given in Table C| The Xc2 meson has a 
natural width of just 2 MeV |2| and is therefore fitted 
with a Gaussian to account for detector resolution. Since 
the XcO has a natural width of 10 MeV [^, comparable 
to the mass resolution (cr « lOMeV/c^), we fit the XcO 
peak with the convolution of Breit-Wigner and Gaussian 
shapes. 

TABLE I: Efficiencies from fits of exclusive MC distributions 
of rrif+p-^ — mp+i- , with statistical uncertainty. 



K*" 0.071 ± 0.001 0.066 ± 0.001 

K*° {KU°) 0.031 ± 0.001 0.020 ± 0.001 

K'^ 0.158 ± 0.001 0.126 ± 0.001 
K*+ {K+Tv°) 0.036 ± 0.001 0.031 ± 0.001 
7^*+ (7^2 7r+) 0.065 ± 0.001 0.062 ± 0.001 

K+ 0.144 ± 0.001 0.117 ± 0.002 



Studies of MC samples show that most of the back- 
ground events in the XcK* channels are due to non- 
resonant (NR) B Xc{ JI4i)Ktt decays. After the NR 
events are removed from the MC background sample, 
the expected background with a genuine Xc ^ J/i^l 
decays is 0.2 ± 0.2 event for the Xc2K*^{K+tt-) and 
Xc2K*+{K+TT°) modes, and 0.0 ± 0.2 for all other chan- 
nels. We correct for the presence of NR decays with the 
following procedure. The TO^+f-..^— to^+^- distribution for 
events in a nearby sideband (1.1 < rriKTr < 1.3 GeV/ c^) 
is subtracted from the distribution for events in the sig- 
nal region (0.85 < < 0.94 GeV/ c^), after scaling the 
sideband distribution by a factor r = 0.26 ± 0.04. The 
quantity r, obtained from MC simulation, is the ratio of 
NR events under the peak to the number in the side- 
band. NR-subtracted distributions of rrii+g-^ — mg+i- 
are shown in Fig.Q] These plots show the presence of the 
factorization-allowed Xci but no significant signals for the 
factorization-suppressed XcO or Xc2- No XcO or Xc2 signal 
is observed in the sideband region. 

TABLE II: Event yields with statistical uncertainties from the 
fits of Fig. El 

Xc2 XcO 

if" {K+TT-) 2.0 ± 1.6 1.7 ± 2.1 

K*° Ik°s-k°) -1.6 ± 4.3 0.5 ± 0.3 

ifg 3.4 ± 1.8 3.9 ± 3.8 

K*+ (A'+tt") -0.5 ± 0.2 1.1 ± 2.2 

K*+ IkI-k+) -1.9 ± 1.2 5.9 ± 3.7 

K+ 3.7 ± 4.4 8.8 ± 6.6 



The branching fractions are computed from BE = 
Ns /{Nsef), where Ns is the number of signal events 
obtained from fitting the to^+^-^ — nii+i- distribution 
C Table lll|) . Nb is the number of produced BB events, 
e is the selection efficiency (Table Pi and / is the prod- 
uct of secondary branching fractions of the B daughters. 
The free parameters in the fits are the size of a constant 
background, the overall scale of mg+g-^ — TO£+£-, and 
the amplitudes of the resonant peaks. The fixed param- 
eters are the XcO natural width, the XcO~Xci and Xc2~Xci 
mass differences (—95.4 and -1-45.7 MeV/c^, respectively) 
all taken from Ref. |^ , and the mass resolution. The mass 
resolution, 10.2 ± 0.4 MeV/ c^, is measured with Xci data 
and is assumed to be the same for the three Xc states. 
Performing such fits to an inclusive T(4S') BB MC 
sample, we verify that the NR events are subtracted cor- 
rectly, and that the proximity of the Xci does not induce 
any significant bias on the measurement of the nearby 

Xc2- 

Based on studies oi B ^ J/ip K* decays [13], the NR 
Kit component appears to be in an S'-wave state, with 
an unknown relative phase (j) with respect to the main 
if* (892) P-wave peak. As no signal is found, the sys- 
tematic uncertainty due to the unknown relative phase 
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mj+,-^ - m,+,- (GeV/c^) 

FIG. 1: Distribution of mi+i~^ — m^+^- for data, with NR 
subtraction for final states of the strange meson (a) K'^n^ , 
(b) T^o^o, (c) Kl (d) T^+ttO, (e) K°7v+, (f) K+ . The fit 
is described in the text. The arrows on plot (f) show the 
expected positions of the XcO and Xc2 peaks. 



lection efficiency therefore depends, to first order, on the 
polarization of the K* population, through the angular 
distribution: 

^ , '^^n = 7 [(1 - cos' Ok') + ^o(3cos2 _ ^ 

(2) 

where Aq is the fraction of longitudinal K* polarization. 
The average efficiency is 

/I (ir 
3 ^eiOK- )d cos Ok' + Ao&, (3) 

where a = | / (I — cos^ 9K')£{dK') sin9K*d9K* , and b = 
J J {3 cos'^ Bk' ~ l)e{dK')sm9K'd9K', where e{9K-') is 
obtained from MC. The values of a and b axe shown in 
Table Iml 

When no signal is observed, as is the case here, the po- 
larization is unknown. We assume an unpolarized decay 
and we estimate the efficiency as (a + 0.56) ± (|&|/vT2). 
The branching fraction measurements reported here are 

TABLE III: Coefficients for the calculation of amplitude- 
dependent average efficiency for the Xa2K* channels (%). 



is estimated here with a MC-based method. The K — tt 
invariant mass is fitted with an amplitude that is the sum 
of a non-relativistic Breit-Wigner and an amplitude with 
a constant phase and the square of which has a quadratic 
dependence on rriKTr- 



rriK* 



rriK-K — ir/2 



b{mK.)e'^ 



(1) 



where a and b are real quantities and mK* = 892 MeV/ . 
The slow variation of the phase of the S wave with rriKT; 
is neglected here. The free parameters in the fit are the 
three degrees of freedom of the quadratic dependence of 
6, the magnitude of the signal, and the relative phase 
(j). As the sideband is dominated by the NR contribu- 
tion, no attempt is made to subtract the few combina- 
torial events. The fact that the phase (j) is unknown is 
dealt with by randomly generating samples of events dis- 
tributed as above for each value of 4>, and applying NR 
subtraction. The number of events N{(j)) thus measured 
is normalized to that obtained with the phase value 4>q 
obtained in the fit. The ratio R = N {(j>) / N {(j)o) shows a 
sinusoidal dependance. The average value is 1.44 with a 
deviation of ±35%, giving an RMS relative uncertainty 
of ±20%, which we will assume as systematic uncertainty 
(due to the interference with the NR component). 

In the case of decays to the tensor Xc2, the efficiency 
depends on the intensity fractions to each of three po- 
larization states. The efficiency is mainly sensitive to 
the value of the K* helicity angle 9k*, because small 
values of 9k' occur for low momentum pions. The se- 



tt b Efficiency 
i^*" {K+n') 8.68 -1.40 7.98± 0.40 
K*" iK°sTv°) 4.25 -1.66 3.43± 0.48 
K*+ {K+n^) 5.05 -1.79 4.16± 0.52 
K*+ {K°siv+) 7.83 -1.84 6.92± 0.53 



affected by the systematic uncertainties described in 
what follows. The relative uncertainty on the number of 
BB events is 1.1%. The secondary branching fractions 
and their uncertainty are taken from Ref. Other es- 
timated uncertainties are: tracking efficiency, 1.3% per 
track added linearly; K° reconstruction, 2.5%; selection 
of the 7 from the Xc decays, 2.5%; 7r° selection, 5.0%; 
PID efficiency, 3.0%. For each mass peak and for AE, 
the uncertainty of the central value and of the width of 
the peaks are measured with the Xci channels. These 
quantities are used to estimate the efficiency uncertainty 
from this source. The ratio of B'^ to B~^ production in 
T{4:S) decays is assumed to be unity. The related uncer- 
tainty is small and is neglected here. A summary of 
the multiplicative contributions to the systematics can 
be found in Table lTvl In addition to these multiplica- 
tive contributions there is a small contribution from the 
uncertainty on r for the NR background subtraction. 

Combining the measurements of the K* sub-modes, 
and with the approximation that the multiplicative ef- 
ficiencies for each K* sub-mode are fully correlated, 
we obtain the branching fractions for the factorization- 
suppressed modes listed in Table Ivl As a cross check, 
the results for the allowed Xci are found to be compati- 
ble with those of a recent analysis |l3| optimized for that 
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TABLE IV: Summary of the multiplicative systematic uncer- 
tainties in percent. The first eight rows are in common to 
decays to XcO and Xc2- 





j\ 7r 




J\ TT 








Number ot _£> s 


i.i 


1 1 
i.i 


1 1 
i.i 


i.i 


i.i 


i.i 


Tracking 


5.2 


2.6 


3.9 


3.9 


3.9 


2.6 


As 








Z.O 




z.o 


1\T dl 1 ^~ T" Q 1 C 


9 ^ 
^ .(J 








9 ^ 


9 ^ 


PID 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


Sample selection 


7.7 


13.1 


11.6 


8.2 


6.5 


6.3 


MC statistics 


1.4 


2.9 


1.7 


1.8 


1.3 


1.3 


S-wave Phase 


20.0 


20.0 


20.0 


20.0 






XcO second. BF 


11.9 


11.9 


11.9 


11.9 


11.9 


11.9 


Total for XcO 


25.4 


28.3 


27.6 


25.5 


14.8 


14.6 


Xc2 second. BF 


8.5 


8.5 


8.5 


8.5 


8.5 


8.5 


Polarization 


5.1 


14.0 


12.4 


7.7 






Total for Xc2 


24.5 


30.5 


29.1 


25.3 


12.2 


12.0 



decay. We obtain upper bounds on the BFs at 90% con- 
fidence level (C.L.) assuming Gaussian statistics for the 
statistical uncertainties and taking into account the sys- 
tematic uncertainties. We have used a Bayesian method 
with uniform prior for positive BF values in the deriva- 
tion of these limits. The upper limits obtained for decays 
to XcO are larger than for Xc2 due to the smaller XcO ra- 
diative BF. For XcqK^ they are compatible with 
the previous measurements ^|^. 

B Xc(o,2)K^*'^ production requires non-factorizable 
contributions. XcoK~^ decays have been previ- 

ously observed. Colangelo et al. [7| explain this with 
rescattering effects and predict a similar rate for B — > 
Xc2K. This is not observed. The upper limits obtained 
for decays to Xc2 are approximately one order of magni- 
tude lower than the branching fractions of the observed 
JS"*" XcoK~^ decays. Furthermore, we find no evidence 
for the decays B Xc(o,2)K* . 
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TABLE V: Upper limits at 90% C.L. and measured branching 
fractions (in pararentheses) in units of 10~^. 



Xc2 X£0 

K*" 0.36 (0.14 ± 0.11 ± 0.14) 7.7 (3.8 ± 2.6 ± 1.5) 

K*+ 0.12 (-0.15 ± 0.05 ± 0.14) 28.6 (13.5 ± 9.6 ± 5.3) 

K+ 0.30 (0.09 ± 0.10 ± 0.11) 8.9 (4.4 ± 3.3 ± 0.7) 

K° 0.41 (0.21 ± 0.11 ± 0.13) 12.4 (5.3 ± 5.0 ± 0.8) 



